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Abstract 



X 



The shell structure underlying shape changes in neutron-rich nuclei near N = 28 has been 



■ investigated by a novel application of the transient field technique to measure the first-excited 

state g factors in '^'^S and ^''S produced as fast radioactive beams. There is a fine balance between 
proton and neutron contributions to the magnetic moments in both nuclei. The g factor of deformed 
^^S does not resemble that of a conventional collective nucleus because spin contributions are more 
important than usual. 

PACS numbers: 21.10.Ky,21.60.Cs,27.30.+t,27.40.+z,25.70.De 
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The fundamental question of how major shell closures, or magic numbers, change in 
neutron-rich nuclei remains unresolved. At present there is conflicting evidence concerning 
the = 28 shell gap in nuclei approaching the neutron dripline. From the measured j3- 
decay halflife it has been suggested that ^48123 is strongly deformed, implying a quenching of 
the N = 28 gap ^J, whereas knockout reactions on ^^Si give evidence for a nearly spherical 
shape P]. Low-excitation level structures and B{E2) values imply that the nearby even 
sulfur isotopes between = 20 and A^ = 28 undergo a transition from spherical at fgS2o, to 
prolate deformed in 16S24 and 16S26, and that the A^ = 28 nucleus 16S28 appears to exhibit 
collectivity of a vibrational character [3,14,15,0]. However the evolution of deformation in 
these nuclei has underlying causes that remain unclear. Some have argued that a weakening 
of the A^ = 28 shell gap is important 0, while others have argued that the effect of adding 
n.„tron,, to the fy. orbit is primarUy to reduce the proton sy,-i„, gap a,rd that a weakenirrg 
of the A^ = 28 shell gap is not needed to explain the observed collectivity near ^"^S [75 . There 
have been several theoretical studies discussing the erosion of the N = 28 shell closure and 
the orrset of defor,nat,o,r (e.g. Re.. BQ and referenees thereirr). 

To resolve questions on the nature and origins of deformation near A^ = 28, we have used 
a novel technique to measure the g factors of the 2i states in fgS22 and ^5824. The g factor, 
or gyromagnetic ratio, is the magnetic moment divided by the angular momentum. The 
existence of deformation in nuclei has long been associated with strong interactions between 
a significant number of valence protons and neutrons, particularly in nuclei near the middle 
of a major shell. Without exception the deformed nuclei studied to date have g factors 
near the hydrodynamical limit, Z/A, reflecting the strong coupling between protons and 
neutrons, and a magnetic moment dominated by the orbital motion of the proton charge 
with small contributions from the intrinsic magnetic moments of either the protons or the 
neutrons. In transitional regions g factors have considerable sensitivity to the proton and/or 
neutron contributions to the state wavef unctions, particularly if the intrinsic spin moments 
of the nucleons come to the fore. 

The present Letter presents the first application of a high-velocity transient-field (HVTF) 
technique to measure the g factors of excited states of fast radioactive beams. In brief. 



intermediate-energy Coulomb excitation is used to excite and align the nuclear states 
of interest. The nucleus is then subjected to the transient field in a higher velocity regime 
than has been used previously for moment measurements, which causes the nuclear spin 
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to precess. Finally, the nuclear precession angle, to which the g factor is proportional, 
is observed via the perturbed 7-ray angular correlation measured using a multi-detector 
array. It is important to note that the transient-field technique has sensitivity to the sign 
of the g factor, which in itself can be a distinguishing characteristic of the proton/neutron 
contributions to the state under study, since the signs of the spin contributions to the proton 
and the neutron g factors are opposite. 

The transient field (TF) is a velocity-dependent magnetic hyperfine interaction expe- 
rienced by the nucleus of a swift ion as it traverses a magnetized ferromagnetic material 
12, Qj- For light ions {Z < 16) traversing iron and gadolinium hosts at high velocity. 



the dependence o 



parametrized 
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_the TF strength on the ion velocity, v, and atomic number, Z, can be 
as 

B,,{v, Z) = AZ''{vlZv^)\-"^'^"l^'"^^\ (1) 



where vq = c/137 is the Bohr velocity. A fit to data for iron hosts yielded A = 1.82(5) T 
with P = 3 10]. The maximum TF strength is reached when the ion velocity matches the 
i^-shell electron velocity, v = Zvq. Since the transient field arises from polarized electrons 
carried by the moving ion its strength falls off as the ion velocity exceeds Zvq and becomes 
fully stripped; a transient-field interaction will not occur for fast radioactive beams with 
energies near 100 MeV/nucleon until most of that energy is removed. 

The experiment was conducted at the Coupled Cyclotron Facility of the National Su- 
perconducting Cyclotron Laboratory at Michigan State University. Secondary beams of ^^S 
and ^°S were produced from 140 MeV/nucleon primary beams directed onto a ~ 1 g/cm^ 
Be fragmentation target at the entrance of the A1900 fragment separator An acrylic 
wedge degrader 971 mg/cm^ thick and a 0.5% momentum slit at the dispersive image of the 
A1900 were employed. The wedge degrader allowed the production of highly pure beams 
and also reduced the secondary beam energy to ~ 40 MeV/nucleon. Further details of the 
radioactive beams are given in Table HI The ^'^S (^°S) measurement ran for 81 (68) hours. 

Figure ^ shows the experimental arrangement. The radioactive beams were delivered 
onto a target which consisted of a 355 mg/cm^ Au layer backed by a 110 mg/cm^ Fe layer 
of dimensions 30 x 30 mm^. The target was held between the pole tips of a compact 
electromagnet that provided a magnetic field of 0.11 T, sufficient to fully magnetize the Fe 
layer. To minimize possible systematic errors, the external magnetic field was automatically 
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TABLE I: Production and properties of radioactive beams. 





Primary beam 






Secondary beam 




Ion 


Intensity 


Ion 


E 


Intensity 


Purity 




(pnA) 




(MeV) 


(pps) 


(%) 


40 


25 


38s 


1547.5 


2 • 10^ 


> 99 


48Ca 


15 


40s 


1582.5 


2 • 10^ 


> 95 




FIG. 1: Schematic view of the experimental arrangement from above. Only the four SeGA detectors 
perpendicular to the magnetic field axis are shown. There are ten other SeGA detectors (see text) 
which are not shown for clarity. 

reversed every 600 s. 

Table |n] summarizes the properties of the 2i states, the key aspects of the energy loss of 
the sulfur beams in the target, the precession results and the extracted g factors. The high- 
Z Au target layer serves to enhance the Coulomb excitation yield and slow the projectiles 
to under 800 MeV, while the thick iron layer results in a long interaction time with the 
transient field, maximizing the spin precession. The sulfur fragments emerge with energies 
in the range from ~ 80 MeV to ~ 200 MeV. Most of this energy spread stems from the 
energy width of the radioactive beam. 

Projectiles scattering forward out of the target were detected with a 15.24 cm diameter 
plastic scintillator phoswich detector placed 79.2 cm downstream of the target position. The 
maximum scattering angle of 5.5° limits the distance of closest approach to near the nuclear 
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interaction radius in both the Au and Fe target layers. Positioning the particle detector 
downstream also lowers the exposure of the 7-ray detectors to the radioactive decay of the 
projectiles. 

To detect de-excitation 7 rays, the target chamber was surrounded by 14 HPGe detectors 
of the Segmented Germanium Array (SeGA) 3|- The SeGA detectors were positioned with 
the crystal centers 24.5 cm from the target position. Six pairs of detectors were fixed at 
symmetric angles (i^, 0) = (29°, 90°), (40°, 131°), (60°, 61°), (139°, 46°), (147°, 143°), and 
(151°, 90°), where 9 is the polar angle with respect to the beam axis and is the azimuthal 
angle measured from the vertical direction, which coincides with the magnetic field axis. 
Each ±6' pair is in a plane that passes through the center of the target; (—6', 0) = (9^ 0+180°). 
Two more detectors were placed at 6* = 90° and 6 = 24°. All 14 detectors were used 
to measure the 7-ray angular correlations concurrently with the precessions. Since the 
precession angles are small, the unperturbed angular correlation can be reconstructed by 
adding the data for the two directions of the applied magnetic field. 

Coincidences between the phoswich particle detector and SeGA were recorded, and 7- 
ray spectra gated on sulfur recoils and corrected for random coincidences were produced. 
Doppler-corrected spectra were also produced using the angular information from the SeGA 
detector segments and the particle energy information from the phoswich detector on an 
event-by-event basis, which is essential because of the spread in particle velocities. Figure |21 
shows examples of the 7-ray spectra. From the measured Doppler shift of the deexcitation 
7 rays in the laboratory frame, the average after-target ion velocities were determined. The 
velocity distribution of the exiting ^''S ions was also measured by shifting the phoswich 
detector by ±15 cm from its normal position and observing the change in the flight times of 
the projectiles. These procedures firmly establish that the sulfur ions were slowed through 
the peak of the TF strength at Zvq into the region where it has been well characterized 

HQ. 

The 2"*" peak areas averaged 925 counts/detector per field direction for ^^S and 400 
counts/detector per field direction for ^°S, in each of the six angle pairs of SeGA detec- 
tors used for extracting the precessions. 



The angular correlation of 7 rays was calculated with the program GKINT [17[ using 
the theory of Coulomb excitation llSJ. Recoil-in- vacuum effects were evaluated based on 

n 

measured charge-state fractions for sulfur ions emerging from iron foils [1^. Good agreement 
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FIG. 2: Doppler-corrected 7-ray spectra at 9 = 40° for ^^S (left) and ^'^S (right). The sulfur and 
iron peaks are labeled. The broad feature to the right of the 1292 keV line is mainly its Doppler 
tail, due to decays within the target. 



1.5 
1.25 

1 ^ 
£ 0.75 

^ 0.5 

0.25 

-180 -90 90 180 -180 -90 90 180 
Qnuc [deg] e„u, [deg] 

FIG. 3: Angular correlations for ^^S and '^'^S in the frame of the projectile nucleus. Data are 
normalized to the calculated angular correlation. The difference in anisotropy stems from the 
alignment produced by Coulomb excitation, which depends on the ratio of E(2"'') to the beam 
velocity. 




was found between the calculated 7-ray angular correlations and the data, as shown in Fig 01 
The experimental nuclear precession angle, A^^exp, was extracted from 7-ray count ratios in 
pairs of detectors at (±6', 0) for both field directions, using standard analysis methods as 
described in |14]. 

An evaluation of A6/g = (— /in/^) / Btfdt is required to extract the g factors. Calcula- 
tions were performed using the code GKINT to take into account the incoming and exiting 
ion velocities, the energy- and angle-dependent Coulomb excitation cross sections in both 
target layers, the excited-state lifetimes, and the parametrization of the TF strength in 
Eq. (^. The results and the g factors extracted are given in Table ITTl These g factor results 
are not very sensitive to the somewhat uncertain behavior of the transient field at the highest 
velocities because (i) the ions spend least time interacting with the TF at high velocity and 
(ii) the TF strength near 2Zvq is very small. Furthermore, the positive g factor in ^^S and 
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TABLE II: Nuclear parameters and reaction kinematics. B{E2)] = B{E2;0^^ — > 2^). (i?j,e) and 
(vi^e) are the ion kinetic energies and velocities at the entrance and exit of the iron layer. The 
effective transient-field interaction time teflf and the spin precession per unit g factor, {A6 / g)cs,ic, 
are evaluated for ions that decay after leaving the target. The experimental g factor is given by 

g = A0cxp/(A0/5)calc. 



Isotope E(2|) 


B{E2) t 


r{2t) 


{E^) {Ee) 


{vi/Zvo) 






(A0/(7)ealc 


A6'cxp g 


(keV) 


(e^fm^) 


(ps) 


(MeV) (MeV) 






(ps) 


(mrad) 


(mrad) 


1292 


235(30) 


4.9 


762 123 


1.75 


0.71 


2.98 


-330 


-43(15) +0.13(5) 


904 


334(36) 


21 


782 145 


1.73 


0.75 


2.99 


-339 


+5(21) -0.02(6) 



the essentially null effect for '^''S are both firm observations, independent of the transient- 
field strength. The experimental uncertainties assigned to the g factors are dominated by 
the statistical errors in the 7-ray count ratios, with a small contribution (10%) from the 
angular correlation added in quadrature. 

Shell model calculations were performed for ^l.'&').n. y^"}-?. and 16^245 and their isotones 
'^gAr2o, i8Ar22 and ^gAr24, using the code OXBASH 20| and the sd-pf model space where 
(for > 20) valence protons are restricted to the sd shell and valence neutrons are restricted 

n 

to the pf shell. The Hamiltonian was that developed in Ref. [21|| for neutron-rich nuclei 
around A^ = 28. These calculations reproduce the energies of the low-excitation states 
to within 200 keV. With standard effective charges of Cp ~ 1.5 and e„ ~ 0.5 they also 
reproduce the measured B{E2) values. The g factors of the 2f states were evaluated using 
the bare nucleon g factors. The calculated g factors are compared with experimental results 
in Fig. 0] Overall the level of agreement between theory and experiment is satisfactory 
given the extreme sensitivity to configuration mixing and the near cancelation of proton 
and neutron contributions in the A^ = 22, 24 isotones (see below). 

In the A^ = 20 isotones, '^^S and ^^Ar, the 2^ state is a pure proton excitation for our 
model space. Two neutrons have been added in the fp shell in the N = 22 isotones ^^S and 
^°Ar. Since ^^S is almost doubly magic, the initial expectation might be that the first-excited 
state of '^'^S would be dominated by the neutron fj/2 configuration weakly coupled to the ^^S 
core, resulting in a g factor near —0.3. In contrast, the near zero theoretical g factor and 
the small but positive experimental g factor require additional proton excitations, which 
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FIG. 4: Theoretical g factors compared with experiment. The previous resuhs for ^^Ar and ^'^Ar 
are from Refs. 

Q and Q, 

respectively. 



indicates strong coupling between protons and neutrons - one of the prerequisites for the 
onset of deformation. For = 22, 24 the shell model predicts a cancelation of the proton and 
neutron contributions to the moment; in terms of g^^ = S'proton + fi'neutron' di^^^) = —0.003 = 
0.298 - 0.301, ^(^°S) = 0.035 = 0.276 - 0.241, c/(^°Ar) = -0.200 = 0.164 - 0.364 and 
5f(^^Ar) = —0.060 = 0.220 — 0.280. The proton contributions to the g factors are dominated 
by the orbital component but the substantial neutron contributions originate entirely with 
the intrinsic spin associated with a dominant occupation of the neutron fi/2 orbit. Some 
tuning of the Hamiltonian may be required to reproduce the sign of the g factor in ^^S, 
which is very sensitive to the separation of the proton si/2 and ^3/2 orbitals, for example. 
The dependence of the g{2'^) in ^°Ar on the basis space, the interaction, and the choice of 
effective nucleon g factors, has been investigated in Ref. [2^. 

As noted above. Coulomb-excitation studies and the level scheme of '^''S suggest that it 
is deformed. Supporting this interpretation, the shell model calculations predict consistent 
intrinsic quadrupole moments when derived from either the B{E2) or the quadrupole mo- 
ment, Q{2f), implying a prolate deformation of /3 ~ -1-0.3, in agreement with the value 
deduced from the experimental B[E2) 0, But the near zero magnetic moment does 
not conform to the usual collective model expectation of (7 ~ Z/A. Since the shell model 
calculations reproduce both the electric and magnetic properties of the 2\ state they give 
insight into the reasons for this unprecedented magnetic behavior in an apparently deformed 
nucleus. The essential difference between the deformed neutron-rich sulfur isotopes and the 
deformed nuclei previously encountered (i.e. either light nuclei with N = Z 01 heavier de- 
formed nuclei) is that the spin contributions to the magnetic moments are relatively more 
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important, especially for the neutrons. In comparison to °S, the 2^ state in the = 20 
nucleus f2Mg2o has a similar excitation energy, hfetime and B{E2). However the g factor 
in ^^Mg might be closer to that of a conventional collective nucleus since the = 20 shell 
closure is known to vanish far from stability. 
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the experiment. This work was supported by NSF grants PHY-01-10253, PHY-99-83810, 
and PHY-02-44453. AES, ANW, and PMD acknowledge travel support from the ANSTO 
AMRF scheme (Austraha). 



[1] S. Grevy, et al, Phys. Lett. B 594, 252 (2004). 
[2] J. Pridmann, et al, Nature 435, 922 (2005). 
[3] H. Scheit, et al, Phys. Rev. Lett. 77, 3967 (1996). 
[4] T. Glasmacher, et al, Phys. Lett. B 395, 163 (1997). 

[5] J. Winger, P. Mantica, R. Ronningen, and M. Caprio, Phys. Rev. C 64, 064318 (2001). 

[6] D. Sohler, et al, Phys. Rev. C 66, 054302 (2002). 

[7] P. D. Cottle and K. W. Kemper, Phys. Rev. C 58, 3761 (1998). 

[8] R. Rodriguez-Guzman, J. Edigo, and L. Robledo, Phys. Rev. C 65, 024304 (2002). 

[9] E. Caurier, F. Nowacki, and A. Poves, Nucl. Phys. A 742, 14 (2004). 
[10] A. E. Stuchbery, Phys. Rev. C 69, 064311 (2004). 
[11] T. Glasmacher, Annu. Rev. Nucl. Sci. 48, 1 (1998). 

[12] N. Benczer-KoUer, M. Hass, and J. Sak, Annu. Rev. Nucl. Sci. 30, 53 (1980). 

[13] K. H. Speidel, O. Kenn, and F. Nowacki, Prog. Part. Nucl. Phys. 49, 91 (2002). 

[14] A. E. Stuchbery, et al, Phys. Lett. B 611, 81 (2005). 

[15] D. J. Morrissey, et al, Nucl. Inst. Meth. Phys. Res. B 204, 90 (2003). 

[16] W. F. Mueller, et al, Nucl. Inst. Meth. Phys. Res. A 466, 492 (2001). 

[17] A. E. Stuchbery, Some notes on the program GKINT: Transient-field g-factor kinematics at 
intermediate energies. Department of Nuclear Physics, The Australian National University, 
report no. ANU-P/1678 (2005). 

[18] C. A. Bertulani, A. E. Stuchbery, T. J. Mertzimekis, and A. D. Davies, Phys. Rev. C 68, 
044609 (2003). 



9 



[19] A. E. Stuchbery, P. M. Davidson, and A. N. Wilson, Nucl. Instr. Meth. Phys. Res. B 243, 
265 (2006). 

[20] B. A. Brown, et al, Oxbash for Windows PC, Michigan State University, report no. MSU- 

NSCL 1289 (2004). 
[21] S. Nummela, et al, Phys. Rev. C 63, 044316 (2001). 
[22] K.-H. Speidel, et al, Phys. Lett. B 632, 207 (2006). 
[23] E. Stefanova, et al, Phys. Rev. C 72, 014309 (2005). 

[24] J. Retamosa, E. Caurier, F. Nowacki, and A. Poves, Phys. Rev. C 55, 1266 (1997). 



10 



